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Abstract. The authors survey recent results in special functions of 
classical analysis and geometric function theory, in particular the circu- 
lar and hyperbolic functions, the gamma function, the elliptic integrals, 
the Gaussian hypergeometric function, power series, and mean values. 



1. Introduction 



The study of quasiconformal maps led the first two authors in their joint 
work with M. K. Vamanamurthy to formulate open problems or questions 
involving special functions |AVVH IAVV3] . During the past two decades, 
many authors have contributed to the solution of these problems. However, 
most of the problems posed in |AVV1| are still open. 

The present paper is the third in a series of surveys by the first two au- 
thors, the previous papers |AVV6| IAVV8] being written jointly with the late 
M. K. Vamanamurthy. The aim of this series of surveys is to review the re- 
sults motivated by the problems in |AVVH IAVV3] and related developments 
during the past two decades. In the first of these we studied classical spe- 
cial functions, and in the next we focused on special functions occurring in 
the distortion theory of quasiconformal maps. Regretfully, Vamanamurthy 
passed away in 2009, and the remaining authors dedicate the present work 
as a tribute to his memory. For an update to the bibliographies of |AVV6| 
and [AVV8J the reader is referred to |AV] , 

In 1993 the following monotone rule was derived |AVV4| Lemma 2.2]. 
Though simple to state and easy to prove by means of the Cauchy Mean 
Value Theorem, this l'Hopital Monotone Rule (LMR) has had wide appli- 
cation to special functions by many authors. Vamanamurthy was especially 
skillful in the application of this rule. We here quote the rule as it was 
restated in |AVV7| Theorem 2]. 

1.1. Theorem. (l'Hopital Monotone Rule). Let -co < a < b < oo, and 
let f,g : [a,b] — > M be continuous functions that are differentiable on (a,b), 
with f(a) = g(a) = or f(b) = g(b) = 0. Assume that g'(x) 7^ for each 
x £ (a,b). If f'/g' is increasing (decreasing) on (a,b), then so is f/g. 



Theorem 11.11 assumes that a and b are finite, but the rule can be extended 
easily by similar methods to the case where a or b is infinite. Clearly, the 
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l'Hopital Monotone Rule is helpful wherever the l'Hopital rules for limits 
are. For instance, Pinelis' note |Pi2| shows the potential of the LMR. As a 
complement to Pinelis' note, the paper |AVV7| contains many applications 
of LMR in calculus. Also the history of LMR is reviewed there. 

In this survey we give an account of the work in the special functions of 
classical analysis and geometric function theory since our second survey. In 
many of these results the LMR was an essential tool. Because of practical 
constraints, we have had to exclude many fine papers and have limited our 
bibliography to those papers most closely connected to our work. 

The aim of our work on special functions has been to solve open problems 
in quasiconformal mapping theory. In particular, we tried to settle Mori's 
conjecture for quasiconformal mappings [ Mo] (see also |LV| p. 68]). For the 
formulation of this problem, let K > 1 be fixed and let M(K) be the least 
constant such that 

\f(zi) - f(z 2 )\ < M(K)\z x ~ z 2 \ 1/K , for all zx,z 2 e B, 

for every If -quasiconformal mapping / : B — )■ B of the unit disk B onto 
itself with /(0) = 0. A. Mori conjectured in 1956 that M(K) < 16 1_1/x . 
This conjecture is still open in 2012. Some of the open problems that we 
found will be discussed in the last section. 

2. Generalizations of Jordan's inequality 

The LMR application list, begun in |AVV7] . led to the Master's thesis of 
M. Visuri, on which |KVV] is based. Furthermore, applications of LMR to 
trigonometric inequalities were given in |KVV] . Numerous further applica- 
tions to trigonometric functions were found by many authors, and some of 
these papers are reviewed in this section and the next. 

By elementary geometric methods one can prove that 

2 sin x 7r 
-< <1, 0<x<-, 

7T X 2 

a result known as Jordan's inequality. In a recent work, R. Klen et al |KVV] 
have obtained the inequalities 



cos^<-<cos^, x€(-y/27/E,y/Z7/E), 

£ X o 

and 

cosh 1//4 x < < cosh 1//2 x, x G (0, 1). 

x 

Inspired by these results, Y.-P. Lv, G.-D. Wang, and Y.-M. Chu |LWC] 
proved that, for a = (log(vr/2))/ log y/2 « 1.30299, 

a in x sm x r, x , . . 

cos 4/3 - < < cos a -, X £ (0,7172), 

2 x 2 

where 4/3 and a are best constants, and that for b = (log sinhl)/(log coshl) « 
0.372168, 

cosh 1 / 3 ^ < Sinh X < cosh b x, x G (0, 1), 

x 

where 1/3 and b are best constants. 

Many authors have generalized or sharpened Jordan's inequality, either by 
replacing the bounds by finite series or hyperbolic functions, or by obtaining 
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analogous results for other functions such as hyperbolic or Bessel functions. 
The comprehensive survey paper by F. Qi et al |QNG| gives a clear picture 
of these developments as of 2009. For example, in 2008 D.-W. Niu et al 
|NHCQ| obtained the sharp inequality 

n n 

- + Y a k (7r 2 - 4* 2 ) fc < — < - + V ^(tt 2 - 4x 2 ) k , 0<x< vr/2, 

k=l k=l 

for each natural number n, with best possible constants and That 
same year S.-H. Wu and H. M. Srivastava |WS2| obtained upper and lower 
estimates on (0, 7r/2] for (sinx)/x that are finite series in powers of (x — 9), 
where 9 G [x,ir/2], while L. Zhu |Z1| obtained bounds as finite series in 
powers of (tt 2 — Ax 2 ). L. Zhu |Z5| obtained bounds for (sinx)/x as finite 
series in powers of (r 2 — x 2 ) for < x < r < tt/2, yielding a new infinite 
series 

oo 

= a n (r 2 - x 2 ) n , for < \x\ < r < vr/2. 

71=0 

S.-J. Yang |Y] showed that a function / admits an infinite series expansions 
of the above type if and only if / is analytic and even. 

In 2011 Z.-H. Huo et al |HNCQ| obtained the following generalization of 
Jordan's inequalities: 



^(e'-xT < — < 22^(6* 



smx smo ,„* ( >t 

1 x) 

X 

k=l fc=l 



for t > 2, n G N, and < x < 9 < tt, where the coefficients and uj\~ are 
defined recursively and are best possible. 

More recently, in 2012, C.-P. Chen and L. Debnath |CD| have proved that, 
for < x < tt/2, 



h(x) < — < f 2 (x), 
x 

where 

= w + ~e~ (7r " {2x) ) + 40^ {7T ~ {2x) ) 

and 

h (x) = - + —n—(* e ~ 2 * " + ^ (n e - (2xf) 2 , 

TT 9 9 

for any 9 > 2, with equality when x = ir/2. 

In a recent work J. Sandor |S3] (see also |S4[ Paper 9]) proved that h(x) = 
[log (x/ sin x)] / log ((sinh a;) /x) is strictly increasing on (0,7r/2). He used this 
result to prove that the best positive constants p and q for which 

sinh x \ p x ( sinh x ^ 

< < 

sin x \ x 

is true are p = 1 and q = [log(7r/2)]/log(sinh(7r/2)/(7r/2)) « 1.18. 

In an unpublished manuscript, C. Barbu and L.-I. Pi§coran [BaPij have 
proved, in particular, that 

(1-x 2 /3)- 1 / 4 <^<1 + ^, XG(0,1). 
x 5 
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M.-K. Kuo [Koj has developed a method of obtaining an increasing se- 
quence of lower bounds and a decreasing sequence of upper bounds for 
(sinx)/x, and he has conjectured that the two sequences converge uniformly 
to {smx)/x. 

Since there is a close connection between the function (sin:c)/cc and the 
Bessel function J\/2{x) (cf. |Z13| ). it is natural for authors to seek analogs of 
the Jordan inequality for Bessel and closely related functions. A. Baricz and 
S.-H. Wu [EarHllEWl] . L. Zhu [ZTH EH] , and D.-W. Niu et al [NCQl have 
produced inequalities of this type. L. Zhu |Z15] has also obtained Jordan- 
type inequalities for {{sm x) / x) p for any p > 0. S.-H. Wu and L. Debnath 
[WDlJ have generalized Jordan's inequality to functions f{x)/x on [0, 9] such 
that / is (n + l)-times differentiable, /(0) = 0, and either n is a positive 
even integer with increasing on [0, 9] or n is a positive odd integer 

with /( n+1 ) decreasing on [0,9]. 



3. Other inequalities involving circular and hyperbolic 

functions 

3.1. Redheffer. In 1968 R. Redheffer [Re] proposed the problem of showing 
that 

, , sirnrx 1 - x 2 

(3.1j > k, tor ail real x, 

TTX 1 + X Z 

or, equivalently, that 

sinx 7r 2 — x 2 

(3.2) > — 7. for all real x. 

x tt z + x 2 

A solution of this problem was provided by J. T. Williams |Ws| . using 
infinite products, who also proved the stronger inequality 

sin7rx 1 — x 2 (1 — x) 2 

TTX 1 + X z X(l + X z ) 

Later, using Erdos-Turan series and harmonic analysis, J.-L. Li and Y.-L. Li 
|LL| proved the double inequality 

(1 — x 2 )(4 — x 2 )(9 — x 2 ) sin7rx 1 — x 2 

— R — A J -\ < < , for < x < 1. 

x e - 2x 4 + 13x 2 + 36 ~ ttx ~ Vl + 3x 4 

They also found a method for obtaining new bounds from old for (sinx)/x, 
but M.-K. Kuo [Ko] gave an example to show that the new bounds are not 
necessarily stronger. In 2003 C.-P. Chen et al |CZQ|, using mathematical 
induction and infinite products, found analogs of the Redheffer inequality 
for cosx: 

IX 2 — 4x 2 7T 

cosx > — 7, — for \x\ < — , 

7r 2 + Ax 2 11 2 

and for hyperbolic functions: 

sinh x tt 2 + x 2 , 7r 2 + 4x 2 . , 7T 

< for < \x \ < 7r; coshx < — ^ — j, tor \x\ < — . 

x ir z — x z 7T 2 — 4x z 2 
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[w+^J ^ C0Sx£ [^TI^) ' ^ x ^ — , 



In 2008, inspired by the inequalities above, L. Zhu and J. -J. Sun [ZSJ 
proved that 

for < x < 

2 

with best possible constants a = 1 and f3 = 7r 2 /16, and 
7r 2 — x 2 \ 7 sin x / 7T 2 — x 2 \ 5 

-o- 2 < < -g- g , for < X < 7T, 

TT Z + X z J X \7T Z + X z J 

with best possible constants 7 = 1 and 5 = 7r 2 /12. They obtained similar 
results for the hyperbolic sine and cosine functions. In 2009 L. Zhu [ZllJ 
showed that 

2 2 \ a ■ / 2 2 \ P 

■k — x \ smx / TT — x z \ 

< < ; , < X < 7T, 



V7T 4 + 3x 4 y x VVvr 4 + 3x 4 , 

holds if and only if a > 7r 2 /6 and /3 < 1, with analogous results for cosx and 
(tanx)/x. In 2009 A. Baricz and S.-H. Wu |EW2j and in 2011 L. Zhu [Z16] 
proved Redheffer-type inequalities for Bessel functions. 

3.2. Cusa-Huygens. The inequality 

sin x 2 + cos x 

< o » < x < tt/2, 

x 3 

was discovered by N. de Cusa in the fifteenth century (cf. |CC1| ). and proved 
rigorously by Huygens |Hs] in the seventeenth century. In 2009 L. Zhu |Z8] 
obtained the following inequalities of Cusa-Huygens type: 

/ smx \ 2 1. 7r 

< 77 + o( COS2; ) ' < x < -, a > 1, 



V x / 3 3 V ' 2 



and 

/ sinh x \ Q 2 1 



. , < - + -(coshx) a , x > 0, a > 1. 

\ x / 3 3 

That same year L. Zhu |Z9| discovered a more general set of inequalities of 
Cusa type, from which many other types of inequalities for circular functions 
can be derived. He proved the following: Let < x < tt/2. If p > 1, then 

(3.3) (l-a) + a(cosx) p < f— ^ < (1 - (3) + /3(cos xf 



if and only if (3 < 1/3 and a > 1 - (2/vr) p . If < p < 4/5, then §3§ holds if 
and only if a > 1/3 and (3 < 1 — (2/ir) p . If p < 0, then the second inequality 
in (|3.3p holds if and only if /? > 1/3. In a later paper [Z 13] L. Zhu obtained 
estimates for (sin x) /x and (sinh x) jx that led to new infinite series for these 
functions. For some similar results see also [WBJ. 

In 2011 C.-P. Chen and W.-S. Cheung [CUT] obtained the sharp Cusa- 
Huygens-type inequality 

2 + cos x \ a sin x / 2 + cos x x 



< < 



for < x < 7r/2, with best possible constants a = (log(7r/2))/log(3/2) 
1.11 and (3 = 1. 
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In 2011 E. Neuman and J. Sandor |NS4| discovered a pair of optimal 
inequalities for hyperbolic and trigonometric functions, proving that, for 
< x < 7r/2, the best positive constants p and q in the inequality 

1 sin x 1 

< < 



(coshx) p x (coshx) 9 

are p = (log(7r/2))/logcosh(-7r/2) ~ 0.49 and q = 1/3, and that for x ^ 
the best positive constants p and q in the inequality 

/ sinhx\ p 2 / sinhx\ 9 

\ X J cos X + 1 \ X J 

are p = 3/2 and q= (log 2)/ log[(sinh(7r/2))/(vr/2)] ra 1.82. 

3.3. Becker-Stark. In 1978 M. Becker and E. L. Stark |BeS] obtained the 
double inequality 

8 tanx tt 2 7r 
< < —r -r, < x < 



TT 2 — Ax 2 X IT 2 — Ax 2 ' 2 

where the constants 8 and tt 2 are best possible. 
In 2008 L. Zhu and J. -J. Sun [ZSJ showed that 

tt 2 + Ax 2 \ ° tan x ( ir 2 + Ax 2 \ 13 



TT 



it 2 -Ax 2 ) - x -\tt 2 -Ax 2 ) ' U<x< 2 ' 

holds if and only if a < tt 2 /2A and j3 > 1. 

In 2010 L. Zhu and J.-K. Hua |ZH| sharpened the Becker-Stark inequality 
by proving that 

■K 2 + ax 2 tanx tt 2 + fix 2 tt 

TT 2 — Ax 2 X TT 2 — Ax 2 ' 2 ' 

where a = 4(8 - tt 2 )/tt 2 -0.76 and /3 = vr 2 /3 - 4« -0.71 are the best 
possible constants. They also developed a systematic method for obtaining 
a sequence of sharp inequalities of this sort. 
In 2011 H.-F. Ge [G] obtained 

tanx 7r 4 1 / 7r 2 

+ l--n < < TT7 o T^7+ 1 



7T 2 - 4X 2 V ^ 2 7 » 12 71-2 - 4 ^ 2 V 12 

for < x < it/2. That same year C.-P. Chen and W.-S. Cheung |CC1] 
proved the sharp Becker-Stark-type inequality 

2 \ a / 1 \ 

tt \ tanx / 7r 

< < 



7T 2 — 4x 2 / x \7r 2 — 4x 2 , 

with best possible constants a = tt 2 /12 ~ 0.82 and ft = 1. 

3.4. Wilker. In 1989 J. Wilker [W] posed the problem of proving that 



/ sinx\ 2 tanx „ tt 

3.4 + > 2, for < x < -, 

\ x J x 2 

and of finding 
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J. Anglesio et al [WSJVAJ showed that the function in H3 . 5 [) is decreasing 
on (0, vr/2), that the value of c in f j 3 . 5 [) is 16/7T 4 , and that, moreover, the 
supremum of the expression in fj3 . 5 [) on (0, vr/2) is 8/45. Hence 

16 o /sinx\ 2 tanx 8 o 

(3.6) 2 + —rx i tanx < [ H <2 + — x d tanx, 

7T 4 ~ \ x J x ~ 45 

for < £ < 7r/2, where 16/7T 4 w 0.164 and 8/45 « 0.178 are best possible 
constants. (Note: |AVV7| erroneously quoted |WSJVA| as saying that the 
function in (|3.5p is increasing.) In 2007 S.-H. Wu and H. M. Srivastava |WS1] 
proved the Wilker-type inequality 

2 



(3.7) 



smx 



+ 



7T 



> 2, for < x < -. 

tan x 2 



However, A. Baricz and J. Sandor |BS| discovered that (|3.7p is implied by 

In 2009 L. Zhu [Z8] generalized (|3.4p and obtained analogs for hyperbolic 
functions, showing that, for < x < 7r/2, a > 1, 

2a 



SIM 



tan x 



and that, for x > 0, a > 1, 



sinhx 



2a 



/ tanh; 



> 



> 



2a 



SIM 



sinh x 



x 



tanx 



>2, 



2o 



+ 



tanhx 



> 2. 



These two results of Zhu are special cases of a recent lemma due to E. 
Neuman [N4, Lemma 2]. 

Using power series, C.-P. Chen and W.-S. Cheung |CC2| obtained the 
following sharper versions of H 3 . 6 [) : 



16 c 

(3.8) x b tan x < 

315 

and 
(3.9) 

104 7 , / sinx 

4725 : 



SIM 



tan x 



2 + 77^ 
45 



< ( — ) x tanx, 

7T 



-x tan x < 



+ 



tanx 



x 



45' 



16 

315" 



< 



x 7 tanx. 



The constants 16/315 « 0.051 and (2/vr) 6 ps 0.067 in ([331) and 104/4725 rj 
0.022 and (2/vr) 8 « 0.027 in (|23|) are best possible. For < x < vr/2, Chen 
and Cheung also obtained upper estimates complementary to 



(— ) : 

V sin x / 



and 



SIM 



SIM 

2 

+ 



H < 2 + — x 6 tan x 

tan x 45 



T— < 2 +77 x ' i + 7T77 x ° tan:E ' 
tan x 45 945 



where the constants 2/45 and 8/945 are best possible. 
In 2012, J. Sandor IS3I has shown that 



sin x sinh x 

Vq > q+l, x / 0, 

x x 



8 G. D. ANDERSON, M. VUORINEN, AND X.-H. ZHANG 

and 

/ sinh x \ q sin x it 

+ > 2, < x < -, 

\ x J x 2 

where q = [log(vr/2)]/log[(sinh(7r/2))/(vr/2)] « 1.18. 

Extensions of the generalized Wilker inequality for Bessel functions were 
obtained by A. Baricz and J. Sandor [BSJ in 2008. 

3.5. Huygens. An older inequality due to C. Huygens [HsJ is similar in form 
to (CUD : 

/o ,^ r, /sinx\ tanx „ , _ . . tt 

(3.10) 2 H > 3, for < \x \ < -, 

\ x J x 2 

and actually implies (J31J (see |NS3] ). In 2009 L. Zhu [Z7j obtained the 
following inequalities of Huygens type: 

, . sin x tan x , . sin x tan x 
(1 -p) hp > 1 > (1 - q) Vq 

XX XX 

for all x S (0, 7r/2) if and only if p > 1/3 and q < 0; 

, sinh x tanh x , , sinh x tanh x 
(1 -p) hp > 1 > (1 - q) h q 

XX XX 

for all x E (0, oo) if and only if p < 1/3 and q > 1; 

1-p +p- > 1 > (1-g + 9 - 

smx tanx smx tanx 

for all x G (0, 7r/2) if and only if p < 1/3 and q > 1 — 2/ir; and 

(ry* ■ y * -'v* • v • 

1 -p)-r-r KPt — 7— > 1 > (1 - V 1- 



sinhx tanhx sinhx tanhx 

for all x E (0, oo) if and only if p > 1/3 and q < 0. 

C.-P. Chen and W.-S. Cheung |CC2| also found sharper versions of (|3.10p . 
as follows: For < x < 7r/2, 

„ 3 o „ /sinx\ tanx „ /2\ 4 , 
(3.11) 3H x 3 tanx<2 M <3+ - x 3 tanx 

20 \ X J X \ 7T / 



and 

,„ 3 c „ / sinx\ tanx 

(3.12) — x 5 tanx<2 



56 V x / x 



o i ^ 4 
3 H x 

20 



, 2\ b , 
< ( — I x tan x, 



where the constants 3/20 = 0.15 and (2/vr) 4 w 0.16 in (l3~Tll and 3/56 w 0.54 
and (2/vr) 6 « 0.067 in (I3TT2]) are best possible. 

Recently Y. Hua |Hu] have proved the following sharp inequalities: For 
< |x| < 7r/2, 

1 q sinx tan(x/2) 80 - 24-7T , 
3 + — x 6 sin x < h 2 -r-^ < 3 H ^ x A sin x, 

40 X x/2 7T 4 

where the constants 1/40 and (80 — 247r)/7r 4 are best possible; and, for i/O, 

3 q sinhx tanhx 3 q 

3 H x tanhx < 2 1 < 3 H x sinhx, 

20 x x 20 

where the constant 3/20 is best possible. 
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3.6. Shafer. The problem of proving 



arctanx > 



3.r 



1 + 2^1 + x 2 '' 



x > 0, 



was proposed by R. E. Shafer in 1966. Solutions were obtained by L. S. 
Grinstein, D. C. B. Marsh, and J. D. E. Konhauser |SGMK| in 1967. In 2011 
C.-P. Chen, W.-S. Cheung, and W.-S. Wang |CCWj found, for each a > 0, 
the largest number b and the smallest number c such that the inequalities 



bx 



1 + ay/l + x 2 



< arctanx < 



cx 



1 + ay/l + x 2 



are valid for all x > 0. Their answer to this question is indicated in the 
following table: 



a 


largest b 


smallest c 


< a < vr/2 


b = -kcl/2 


c = 1 + a 


7r/2 < a < 2/(tt - 2) 


b = 4(a 2 - I) /a 2 


c = 1 + a 


2/(tt - 2) < a < 2 


b = 4(a 2 - I) /a 2 


c = ira/2 


2 < a < oo 


b = l + a 


c = ira/2 



In 1974, in a numerical analytical context |jSh2j, R. E. Shafer presented 
the inequality 



arctanrc > 



3 + ^25 + (80/3)x 2 



x > 0, 



which he later proved analytically |Sh3j . In [Z2J L. Zhu proved that the 
constant 80/3 in Shafer 's inequality is best possible, and also obtained the 
complementary inequality 



arctanx < 



8.r 



3 + x/25 + (256/tt 2 



x > 0, 



where 256/7T 2 is the best possible constant. 

3.7. Fink. In |M7| p. 247], there is a lower bound for arcsinx on [0,1] 
that is similar to Shafer's for arctanx. In 1995 A. M. Fink |F] supplied a 
complementary upper bound. The resulting double inequality is 



(3.13) 



3x 



2 + x/l - x 2 



< arcsinx < 



7TX 



2 + VT 



< x < 1, 



x^ 



and both numerator constants are best possible. Further refinements on 
these inequalities, along with analogous ones for arcsinhx were obtained by 
L. Zhu [ZH] and by W.-H. Pan with L. Zhu [P ZJ . We note that, for < x < 1, 
this double inequality is equivalent to 

2 + cos x sin x 2 + cos x ir 
< < ^ , 0<x<-, 

7T X 3 2 



in which the second relation is the Cusa inequality. 



10 



G. D. ANDERSON, M. VUORINEN, AND X.-H. ZHANG 



3.8. Carlson. In 1970 B. C. Carlson |C3j (1.14)] proved the inequality 
(3.14) 



6VT 



< arccosx < —t-_ , , ,„ , < x < 1. 



2yj2 + y/TTx~ " " (l + x)V6 

In 2012, seeking to sharpen and generalize (|3.14p . C.-P. Chen and C. Mortici 
[CM2J determined, for each fixed c > 0, the largest number a and smallest 
number b such that the double inequality 



ay/1 — x b\/l — x 

< arccosx < 



C + \/l + X C + \/l + X 

is valid for all x S [0, 1]. Their answer to this question is indicated in the 
following table: 



c 


largest a 


smallest 6 


< x < (2vr - 4)/(4 - vr) 


(l + o)7r/2 


2 + ^/2a 


(2vr - 4)/(4 - vr) < x < (4 - ir)/(ir - 2y/2) 


8(a 2 -2)/a 2 


2 + ^2a 


(4 - vr)/(vr - 2y/2) < x < 2^2 


4(a 2 -l)/a 2 


(l + o)7r/2 


2\/2 < x < oo 


2 + ^2a 


(1 + a)7r/2 



These authors also proved that, for all x £ [0, 1], the inequalities 

< arccosx < 



a + (l + x)V6- - 6+(i + aj )i/6' 

hold on [0, 1], with best constants a = (2\/4 — tt)/t: ~ 0.01 and 6 = 0. 

Moreover, in view of the right side of (|3.14p . in 2011 C.-P. Chen, W.-S. 
Cheung, and W.-S. Wang [CCWJ considered functions of the form 



r(l 



{l + x)i 

on [0, 1], and determined the values of p, q, r such that f{x) is the best 3rd- 
order approximation of arccos x in a neighborhood of the origin. The answer 
is that, for p = (it + 2) / 7r 2 , q = (it — 2) /vr 2 , r = 7r/2, one has 



lim 



arccos x — f(x) tt 2 



X 3 67T 2 

With the values of p, q, r stated above, the authors were led to a new lower 
bound for arccos: 

(vr/2)(l - x)(* +2 y* 2 
arccosx > —, — ^77-^5 — , < x < 1. 

(l + x)^- 2 )/^ 2 

3.9. Lazarevic. In |L] I. Lazarevic proved that, for 1 / 0, 

/ sinh x\ 9 



\ x 



> cosh x 



if and only if q > 3. L. Zhu improved upon this inequality in jZlOj . by 
showing that ifp>lorp<8/15 then 

smh x \ 9 

> p + (1 — p) coshx 



for all x > if and only if q > 3(1 — p). For some similar results see also 
|WBj . 
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In 2008 A. Baricz |Bar7| extended the Lazarevic inequality to modified 
Bessel functions and also deduced some Turan- and Lazarevic-type inequal- 
ities for the confluent hypergeometric functions. 



3.10. Neuman. E. Neuman [N5J has recently established several inequal- 
ities involving new combinations of circular and hyperbolic functions. In 
particular, he has proved that if x ^ 0, then 

. 9 /q sinhx l + 2coshx 

coshx 2/3 < . < , 

arcsm(tann x) 3 

,i/9 o .i/q sinhx (cosh 2a;) 1 / 2 + 2 cosh x 

(cosh 2x 1/2 cosh 2 x 1/3 < — — < i '- , 

arcsmh(tanhx) 3 

and 

.i/q sinhx 2(cosh 2x) 1//2 + coshx 

cosh 2x cosh x] 1/3 < ; — < -i '- . 

arctan(tanhx) 3 

4. EULER'S GAMMA FUNCTION 
For > the gamma function is defined by 



|>00 

r(*)= / t z - l e- l dt, 

Jo 



and the definition is extended by analytic continuation to the entire complex 
plane minus the set of nonpositive integers. This function, discovered by 
Leonhard Euler in 1729, is a natural generalization of the factorial, because 
of the functional identity 

r(* + i) = zT{z). 

The gamma function is one of the best-known and most important special 
functions in mathematics, and has been studied intensively. 

We begin our treatment of this subject by considering an important special 
constant discovered by Euler and related to the gamma function. 

4.1. The Euler-Mascheroni constant and harmonic numbers. The 

Euler- M as cheroni constant 7 = 0.5772156649... is defined as 

(4.1) 7 = lim 7„, 

n— i>oo 

where j n = H n — log n, n E N, and where H n are the harmonic numbers 

l-x r > 



n 1 r 1 1 



-dx. 



x 



k=l 

The number 7 is one of the most important constants in mathematics, and 
is useful in analysis, probability theory, number theory, and other branches 
of pure and applied mathematics. The numerical value of 7 is known to 
29, 844, 489, 545 decimal places, thanks to computation by Yee and Chan in 
2009 |Ye] (see [CMl p. 273]). 

The sequence 7„ converges very slowly to 7, namely with order 1/n. By 
replacing logra in this sequence by log(n+l/2), D. W. DeTemple jD] obtained 
quadratic convergence (see also |Cnl| ). In |M4| C. Mortici made a careful 
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study of how convergence is affected by changes in the logarithm term. He 
introduced new sequences 

Pin) 
Q[n) 

where P and Q are polynomials with leading coefficient 1 and deg P— deg Q = 
1. By judicious choice of the degrees and coefficients of P and Q he was able 
to produce sequences M n tending to 7 with convergence of order 1/n and 
1/ra 6 . He also gave a recipe for obtaining sequences converging to 7 with 
order l/n 2k+2 , where k is any positive integer. This study is based on the 
author's lemma, proved in |M5| . that connects the rate of convergence of a 
convergent sequence {x n } to that of the sequence {x n — x n+ i}. 

In 1997 T. Negoi [M] showed that if T n = H n - log(n + 1/2 + l/(24n)), 
then T n + [4n 3 ] -1 is strictly decreasing to 7 and T n + [48(n + l)] -3 is strictly 
increasing to 7, so that [48(n + l)]~ 3 < 7~T n < [48n 3 ] _1 . In 2011 C.-P. Chen 
[Cn2| established sharper bounds for 7 — T n , by using a lemma of Mortici 
[M5] , 

Using another approach, in 2011 E. Chlebus [Chi] developed a recursive 
scheme for modifying the sequence H n — log n to accelerate the convergence 
to 7 to any desired order. The first step in Chlebus' scheme is equivalent to 
the DeTemple [D] approximation, while the next step yields a sequence that 
closely resembles the one due to T. Negoi |Ni| . 

In |A17| H. Alzer studied the harmonic numbers (|4.2p . obtaining several 
new inequalities for them. In particular, for n > 2, he proved that 

logflog n + 7) 1/n i/(n+i) fl log(logn + 7 ) 

^2 - H n ~ H n+l < P ^2 ' 

where a = (6\/6 - 2 \/396)/ (3 log (log 2 + 7)) w 0.014 and (3 = 1 are the best 
possible constants and 7 is the Euler-Mascheroni constant. 

4.2. Estimates for the gamma function. In |AVVH Lemma 2.39] G. D. 
Anderson, M. Vamanamurthy, and M. Vuorinen proved that 

(4.3, Ita 1 ° gr ,(i + 1 ) = I 

x-5-oo x log X 2 

and that the function (logr(l + x/2))/x is strictly increasing from [2, 00) 
onto [0,oo). In |AQ| G. D. Anderson and S.-L. Qiu showed that the function 
(logr(x + l))/(xlogx) is strictly increasing from (l,oo) onto (1—7,1), where 
7 is the Euler-Mascheroni constant defined by (|4.ip . thereby obtaining the 
strict inequalities 

(4.4) x (i-7)*-i < r(x ) < x x-l : x > 1. 

They also conjectured that the function (logT(x + l))/(xlogx) is concave 
on (l,oo), and this conjecture was proved by A. Elbert and A. Laforgia in 
[ELl Section 3]. 

Later H. Alzer |A12| was able to extend (|4.4p by proving that, for x G (0, 1), 

(4.5) x a(*-l)-7 < T(x) < X^- 1 ^, 

with best possible constants a = 1 — 7 = 0.42278 . . . and /3 = (n 2 /6 — 7)/2 = 
0.53385 .... For x E (l,oo) H. Alzer was able to sharpen (|4.4p by showing 
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that (|4,5p holds with best possible constants a = (vr 2 /6 — r y)/2 ~ 0.534 
and (3 = 1. His principal new tool was the convolution theorem for Laplace 
transforms. 

Another type of approximation for T(x) was derived by P. Ivady jl] in 
2009: 

, , x 2 + 1 . . x 2 + 2 

(4.6 -<T(x + l < -, 0<x<l. 

v ' x+l y ' x + 2 

In 2011 J.-L. Zhao, B.-N. Guo, and F. Qi |ZGQ| simplified and sharpened 
,6p by proving that the function 

Q(x) ^T(x + 1) 



log(x 2 + 1) — log (a; + 1) 

is strictly increasing from (0,1) onto (7,2(1 — 7)), where 7 is the Euler- 
Mascheroni constant. As a consequence, they proved that 

+ ' ^ <T(x + l)<( ^±lY , < x < 1, 



x+l J \x+l 
with best possible constants a = 2(1 — 7) and /3 = 7. 



Very recently C. Mortici |M8| has determined by numerical experiments 
that the upper estimate in (|4,4p is a better approximation for T(x) than 
the lower one when x is very large. Hence, he has sought estimates of the 
form T(x) ~ x a ( x \ where a(x) is close to x — 1 as x approaches infinity. For 
example, he proves that 

x (x-l)a(x) < T ^ < x {x-l)b(x)^ x > 20, 

where a(x) = 1 — 1/logx + l/(2x) — (1 — (log 27r)/2)/(xlog x) and b(x) = 
1 — 1/logx + l/(2x). The left inequality is valid for x > 2. Mortici has 
also obtained a pair of sharper inequalities of this type, valid for x > 2, and 
has showed how lower and upper estimates of any desired accuracy may be 
obtained. His proofs are based on an approximation for logT(x) in terms 
of series involving Bernoulli numbers |AAR| p. 29] and on truncations of 
an asymptotic series for the function (log T(x))/ ((x — 1) log x). These results 
provide improvements of (|4.4[) . 

4.3. Factorials and Stirling's formula. The well-known Stirling's for- 
mula for n\, 

-J V2^i, 

discovered by the precocious home-schooled and largely self-taught eigh- 
teenth century Scottish mathematician James Stirling, approximates n\ asymp- 
totically in the sense that 

n! 

lim — = 1. 

Because of the importance of this formula in probability and statistics, num- 
ber theory, and scientific computations, several authors have sought to re- 
place (|4.7p by a simple sequence that approximates n! more closely (see the 
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discussions in |Bat3] and |Bat4| ). For example, W. Burnside [BuJ proved in 
1917 that 

'n + l/2^ +1/2 



(4.8) n\ ~ fa = V2tt 
that is, lim (n!//3 n ) = 1. In 2008 N. Batir |Bat3] determined that the best 

n— >oo 

constants a and b such that 

(4.9) = — < n! < — 

yn — a yn — b 

are a = 1 — 2-7re -2 ~ 0.1497 and b = 1/6 ~ 0.1667. Batir offers a numerical 
table illustrating that his upper bound formula n n+1 e~ n a/27t / yn — 1/6 gives 
much better approximations to n! than does either (|4.7p or (|4.8p . 

In a later paper |Bat4| N. Batir observed that many of the improvements 
of Stirling's formula take the form 

(4.10) n \„ e -*fa±*\ ^2n{n + b) 

for some real numbers a and b. Batir sought the pair of constants a and 
b that would make f|4. 10[> optimal. He proved that the best pairs (a, b) are 
(a\,bi) and (02,62); where 

ai = - + - - -a/6- A 2 + 4/A « 0.54032, 61 = a? + 1/6 0.45861 
3 6 6 

and 

A 1 



«2 = 77 + 77 + -a/6- A 2 + 4/A w 0.95011, 6 2 = ai + 1/6 « 1-06937, 
3 o o 



where A = a/2 + 2 2 /3 + 2^/3 ~ 2.47128 and a x and 02 are the real roots of 
the quartic equation 3x — 4x 3 + x 2 + 1/12 = 0. 

S. Ramanujan [Raj sought to improve Stirling's formula (|4.7p by replacing 
-/2n in the formula by the sixth root of a cubic polynomial in n: 



( n \ n « / I 
(4.11) r(n + l)« v / vr - \ 8n 3 + An 2 + n + —. 



. e / V 30 
In this connection there appears in the record also his double inequality, for 
x > 1, 



(4.12) ^8x3 + 4^ + x + _L < ^ + *] < y 8 x3 + 4x 2 + x + ^. 

Motivated by this record, the authors of |AVV5| defined the function h(x) = 
u(x) e — (8x 3 + 4x 2 + x), where u(x) = (e/x) x T(x + 1) j ^pn., and conjectured 
that h(x) is increasing from (l,oo) into (1/100, 1/30). In 2001 E. A. Karat- 
suba [KaJ settled this conjecture by showing that h(x) is increasing from 
[1, 00) onto [h(l), 1/30), where h(l) = e 6 /ir 3 - 13 » 0.011. 

In an unpublished document, E. A. Karatsuba suggested modifying Ra- 
manujan's approximation formula (|4.1ip by replacing the radical with the 
2fc th root of a polynomial of degree k, and determining the best such asymp- 
totic approximation. Such a program was partially realized by C. Mortici 
|M7| in 2011, who proposed formula (|4.13p below for k = 4, but the more 
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general problem suggested by Karatsuba remains an open problem. Mortici's 
proposed Ramanujan-type asymptotic approximation is as follows: 



n\n J 32 32 176 

7r | — \ Ion 4 H n 6 H n z H n 

e/ V 3 9 405 



(4.13) r(n + l)f 
In connection with (|4.13j) . lie denned the function 



128 
1215 



g(x) = u(xf 



32 



32 



16x 4 + —x* + —x z + 

3 » 



176 

405 : 



where u(x) = (e/x) x T(x + and proved that g{x) is strictly decreas- 

ing from [3, oo) onto (g(oo), <?(3)], where 5(00) = —128/1215 « —0.105 and 
g(3) = 256e 24 /(43046721vr 4 ) - 218336/135 » -0.088. Mortici's method for 
proving monotonicity was much simpler than Karatsuba's, because he em- 
ployed an excellent result of H. Alzer |A11| concerning complete monotonicity 
(see section l4~6l below for definitions). Mortici claimed that his method would 
also simplify Karatsuba's proof in |Ka| . Finally, he proved that, for x > 3, 

R(x,a)< F( £~tll <R(x,(3), 

7T - 



where R(x,t) = ^/16x 4 + fx 3 + fx 2 + Jggz - t, and a = 128/1215, f3 = 
g(3) are the best possible constants. 

In 2012 M. Mahmoud, M. A. Alghamdi, and R. P. Agarwal |MAA] deduced 
a new family of upper bounds for T(n+ 1) of the form 



r(n + 1) < V2-KV [ - 



n G N, 



Ml 



Tfi] 

n 



1 



2m + 3 



1 m 

j- + E 

An ^ 



k=l 



2m -2k + 2, 
2k + 1 



2~ 2fc C(2A;,n + 1/2) 



n G N, 



where ^ is the Hurwitz zeta function 



k=0 



(k + q y 



These upper bounds improve Mortici's inequality (|4.13p . 

4.4. Volume of the unit ball. The volume Q n of the unit ball in 
given in terms of the gamma function by the formula 

W 2 



IS 



r(ra/2 + l)' 



n G N. 



Whereas the volume of the unit cube is 1 in all dimensions, the numbers 
Q n strictly increase to the maximum = 87r 2 /15 and then strictly decrease 
to as n — > 00 (cf. |BH( p. 264]). G. D. Anderson, M. K. Vamanamurthy, 
and M. Vuorinen |AVV1| proved that Qn is strictly decreasing, and that 
the series Y2n=2 g ^ s convergent. In |AQ| G. D. Anderson and S.-L. 
Qiu proved that Q,l/^ nlogn ^ is strictly decreasing with limit e -1 / 2 as n — > 00. 
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In 2008 H. Alzer published a collection of new inequalities for combinations 
of different dimensions and powers of £l n [A16, Section 3]. We quote several 
of them below: 

(2ire) n / 2 (2ne) n / 2 

(4.14) «^^<(- + l)^-nfi n+1 <6^ w , n>l, 

where the best possible constants are a = (4-97r/8)(2/( 7 re)) 1 / 2 /e = 0.0829 . . . 
and b = ir~ 1 / 2 = 0.5641...; 

(4.15) o^^-IV^Kfi^, n>2, 

with best possible constant factors a = (4/e 2 )(l — 8/(37r)) = 0.0818 . . . and 
6 = 1/(2*-)= 0.1591...; 

(4-16) ° < < * n>2, 

y/n S2 n _i + S2 n+ i y/n 

with best possible constants a = 3\/2~7r / (6 + Air) = 0.7178 . . . and b = \/2tt = 
2.5066...; and 

(4.17) — <(n + l) — n- < n>2, 

y/n il n iZ n _i y/n 



with best possible constants a = (4 — 7r)\/2 = 1.2139 . . . and 6 = y[2/n )2 = 
1.2533.... 

Alzer's work in [A16J includes a number of new results about the gamma 
function and its derivatives. 

In 2010 C. Mortici [M2], improving on some earlier work of H. Alzer |A13| 
Theorem 1], obtained, for n > 1 on the left and for n > 4 on the right, 



2vr " n n/ ±? +1) V2T 



"n+l 



where a = 64 • 720 11 / 12 • 2 1 / 22 /(103 95 • vr 5 / 11 ) = 1.5714. . . . He sharpened 
work of H. Alzer |A13| Theorem 2] and S.-L. Qiu and M. Vuorinen |QV3| in 
the following result, valid for n > 1: 



2n + l fi„_i 2n + l 1 
<^^<\—. + 



4tt Q n V 47r 167rn 

C. Mortici also proved, in |M2| Theorem 4], that, for n > 4, 

i\^~^ n 2 ( 1\3 

This result improves a similar one by H. Alzer |A13| Theorem 3, valid for 
n > 1], where the exponent on the left is the constant 2 — log 2 7r. 



4.5. Digamma and polygamma functions. The logarithmic derivative 
of the gamma function, ip(x) = ^logT(x) = r'(x)/r(a;), is known as the 
digamma function. Its derivatives ip( n \n > 1, are known as the polygamma 
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functions ip n . These functions have the following representations [AbSj pp. 
258-260] for x > and each natural number n: 

r-oo e -t _ e -xt x 



' n(x + n) 

n=l ' 



and 

poo 4.n„—xt 



poo in,.- -xt. 

v„(x) = (-ir +i / -— ^ = (-ir +i n! y> + /c)-"- 1 . 

Several researchers have studied the properties of these functions. In 2007, 
refining the left inequality in |A14| Theorem 4.8], N. Batir |Batf] obtained 
estimates for ip n in terms of i/j or tp^, with k < n. In particular, he proved, 
for x > and n G N: 

(n — 1)! exp (—ntp(x + 1/2)) < |^n(x)| < (n — 1)! exp (— mp(x)) , 

and, for 1 < k < n — 1, x > 0, 



He also proved, for example, the difference formula 

< ((-ir~Vn(x + i)) _1/n - ((-ir-Vn(x))~ 1/n < /?, 



where a = (rt!£(n + 1)) x / n and /3 = ((n — 1)!) 1 l n are best possible, and 
the sharp estimates 

-7 < V(x) + log(e 1/x - 1) < 0, 

where 7 is the Euler-Mascheroni constant. 

In 2010 C. Mortici [M3J proved the following estimates, for x > and 
n > 1, refining work of B.-N. Guo, C.-P. Chen, and F. Qi \CKXj\: 



1 (n + 3)! 

< IV>n(V)| 



720 x n + 4 



(n-1)! 1 n! ; 1 (n + 1)! 

sc" h 2rr n+1 + 12 x n + 2 



< 0. 



4.6. Completely monotonic functions. A function / is said to be com- 
pletely monotonic on an interval I if (-l) n /( n )(a:) > for all a; € J and all 
nonnegative integers n. If this inequality is strict, then / is called strictly 
completely monotonic. Such functions occur in probability theory, numerical 
analysis, and other areas. Some of the most important completely mono- 
tonic functions are the gamma function and the digamma and polygamma 
functions. The Hausdorff-Bernstein-Widder theorem [Wi, Theorem 12b, p. 
161] states that / is completely monotonic on [0, 00) if and only if there is a 
non-negative measure [i on [0, 00) such that 



00 

xt , 



f(x) = / e- xt d»(t) 
Jo 

for all x > 0. There is a well-written introduction to completely monotonic 
functions in [MSJ. 
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In 2008 N. Batir [Bat2] proved that the following function F a (x) related to 
the gamma function is completely monotonic on (0, oo) if and only if a > 1/4 
and that — F a (x) is completely monotonic if and only if a < 0: 

111 

F a (x) = log - x log x + x - - log(27r) + -ip[x) + 



2 m ' 2 rw 6(x-a)' 
As a corollary he was able to prove, for x > 0, the inequality 

( 1 ,t \ 1 "\ r ( x ) ( 1 ,i s 1 

exp —-ip{x) — — r < -j= < exp —-tp(x) 



2 rv ' 6(x-a)J x x e~ x V2^ \ T y ' 6(x-/3). 

with best constants a = 1/4 and (5 = 0, improving his earlier work with H. 
Alzer PPE] , 

In 2010 C. Mortici [M3J showed that for every n > 1, the functions /, g : 
(0, oo) — > R given by 

f(x) = W (x)|- (re - 1)! -l^-l (n + 1)! +J- (n + 3)! 

/ W _| V nW| %n 2x „ +1 12 xn+2 + 72Q xn+4 

and 

, . (n-1)! In! 1 (n + 1)! . . , .. 

are completely monotonic on (0, oo). As a corollary, since f(x) and g{x) are 
positive, he obtained estimates for |^> n (x)| as finite series in negative powers 
of x. 

G. D. Anderson and S.-L. Qiu |AQ|, as well as some other authors (see 
[AdAlJ), have studied the monotonicity properties of the function f(x) = 
(logr(x + l)) fx. In 2011 J. A. Adell and H. Alzer |AdAl| proved that /' is 
completely monotonic on (— l,oo). 

In the course of pursuing research inspired by |AQ| and |AVV1| (see 
[BP4J), in 2012 H. Alzer [A16J discussed properties of the function 

f(x) = ( 1 - l °f X r ) xlogx, 
V log(l+x)y 

which F. Qi and B.-N. Guo |QG| later conjectured to be completely mono- 
tonic on (0, oo). In |BP4| C. Berg and H. L. Pedersen proved this conjecture. 

In 2001 C. Berg and H. L. Pedersen |BPlj proved that the derivative of 
the function 

xiogx 

is completely monotonic (see also [BP2]). This result extends work of |AQ| 
and [EL]. Very recently, C. Berg and H. L. Pedersen [BP3J showed that the 
function 

xlog(ax) 

is a Pick function when a > 1, that is, it extends to a holomorphic function 
mapping the upper half plane into itself. The authors also considered the 
function 

/ n x/2 \ 1/0" tog*) 

= ^(l + x/2) J 

and proved that log/(x + 1) is a Stieltjes function and hence that f(x + 1) 
is completely monotonic on (0, oo). 
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5. The hypergeometric function and Elliptic integrals 
The classical hypergeometric function is defined by 

1-1/ 1-. / (a,n)(b,n) x n 

F{a,b;c;x) = 2 Fi{a,b;c;x) = > — -, \x\ < 1, 

(c,n) n! 

where (a, n) = a(a+l)(a+2) • • • (a+n— 1) for n G N, and (a, 0) = 1 for a ^ 0. 
This function is so general that for proper choice of the parameters a, b, c, one 
obtains logarithms, trigonometric functions, inverse trigonometric functions, 
elliptic integrals, or polynomials of Chebyshev, Legendre, Gegenbauer, or 
Jacobi, and so on (see |AbS| ch. 15]). 



5.1. Hypergeometric functions. The Bernoulli inequality [Ml] P- 34] may 

be written as 

(5.1) log(l + ct) < clog(l + t), 

where c > 1, t > 0. In [KMVJ some Bernoulli- type inequalities have been 
obtained. 

It is well known that in the zero-balanced case c = a+6 the hypergeometric 
function F(a, b; c; x) has a logarithmic singularity at x = 1 (cf. |AVV5| 
Theorem 1.19(6)]). Moreover, as a special case |AbS| 15.1.3], 

(5.2) xF(l,l,2;x) = log 



1 — x 

Because of this connection, M. Vuorinen and his collaborators [KMSVJ 
have generalized versions of (|5.ip to a wide class of hypergeometric func- 
tions. In the course of their investigation they have studied monotonicity and 
convexity /concavity properties of such functions. For example, for positive 
a, b let g(x) = xF(a, b;a + b;x), x G (0, 1). These authors have proved that 
G{x) = \ogg{e x / (l + e x )) is concave on (—00,00) if and only if l/a + 1/6 > 1. 
And they have shown that, for fixed a,b G (0, 1] and for x G (0, 1), p > 0, 
the function 

F ( a,b;a + b; 



1 + xP V 1 + x'P 

is increasing in p. In particular, 

F (a,b;a + b; v ' - ) < ( — — F [a,b;a + b; 



1 + v 7 ? V 'l + v^/ V 1 + r V'' '! + r 

Motivated by the asymptotic behavior of F(:r) = F(a,b;c;x) as x — > 1 , 
S. Simic and M. Vuorinen have carried the above work further in jSVlj . 
finding best possible bounds, when a,b,c > and < x,y < 1, for the 
quotient and difference 

F(a,b;c;x) + F(a,b;c;y) 

— — — ■ — , F(x )+F y) -F {x + y-xy . 

F(a,b; c;x + y-xy) 
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5.2. Complete elliptic integrals. For < r < 1, the complete elliptic 
integrals of the first and second kind are defined as 

n / 2 dt f 1 dt 



V'l-rWt Jo x /(l-t2)(i_ r 2 t 2 ) 



(5.3) X(r) = 
and 

r /2 i f 1 1 1 - r 2 t 2 

(5.4) £(r) = / VI - r 2 sin 2 t = / \ - =-dt, 



1-t 2 

respectively. Letting r' = y/l — r 2 , we often denote 

OC'(r) = X(r'), £'(r) = £(/). 
These elliptic integrals have the hypergeometric series representations 

(5.5) 3C(r) = |F(i,I;l;r 2 ), £ = ~F (I, -±; 1; r 2 ) . 

5.3. The Landen identities. The functions % and £ satisfy the following 
identities due to Landen (BFj 163.01, 164.02]: 

% (tt^) = (1 + r) x(r) ' 30 (i+t) = + r) 3C ' (r) ' 

/ 2^? \ _ 2£(r) -r /2 3C(r) /l - r\ _ £'(r) + r3C'(r) 
Vl + ry~ 1 + r ' V 1 + r / _ ! + r 

Using Landen 's transformation formulas, we have the following identities 
[VZl Lemma 2.8]: For r G (0, 1), let t = (1 - r)/(l + r). Then 

3C(t 2 ) = (1+r)2 ^'(r 2 ), X'(£ 2 ) = (1 + r) 2 X{r 2 ), 



lit 2 ) 



4 

£'(r 2 ) + (r + r 2 + r 3 )0C'(r 2 ) 



1 + r 



,2 



2 = 4£(r 2 )-(3-2r 2 -r 4 )3C(r 2 ) 
{ ' (1 + r) 2 

Generalizing a Landen identity, S. Simic and M. Vuorinen |SV2| have 
determined the precise regions in the aft-plane for which a Landen inequality 
holds for zero-balanced hypergeometric functions. They proved that for all 
a, b > with ab < 1/4 the inequality 

F (a, 6; a + 6; - ( X + r ) F ( a > 6 ' a + b ' r<1 ) 

holds for r E (0,1), while for a, 6 > with 1/a + 1/6 < 4, the following 
reversed inequality is true for each r £ (0, 1) : 

F (a, 6; a + 6; ^ j > (1 + r)F (a, 6; a + 6; r 2 ) . 

In the rest of the a6-plane neither of these inequalities holds for all r £ (0, 1). 
These authors have also obtained sharp bounds for the quotient 

(l + r)F(a,6;a + 6; r 2 ) 
F(a,6;a + 6;4r/(l + r) 2 ) 
in certain regions of the a6-plane. 
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Some earlier results on Landen inequalities for hypergeometric functions 
can be found in [QVl]. Recently, A. B aricz obtained Landen-type inequali- 
ties for generalized Bessel functions [Barl, Bar9j. 

5.4. Legendre's relation and generalizations. It is well known that the 
complete elliptic integrals satisfy the Legendre relation |BF| 110.10]: 

CJV -rCJv — JVJv — — . 

2 

This relation has been generalized in various ways. E. B. Elliott [E] proved 
the identity 

T(l + A + /i)r(l + p. + v) 



-F1-F2 + -P3-F4 — -^bi 7 ^ 
where 



r(A + /i + z^ + |)r(i + 



F x = F{\ + \-\-v;l + \ + fi-x), F 2 = F(±-A, i + i/;l + /i + i/;l-x), 

^3 =F(| + A, i-z/;l + A + /i;x), F 4 = - A, § + z/; 1 + + ^; 1 - zr). 

Elliott proved this formula by a clever change of variables in multiple inte- 
grals. Another proof, based on properties of the hypergeometric differential 
equation, was suggested without details in |AAR( p. 138], and the missing 
details were provided in [AVV6J. It is easy to see that Elliott's formula 
reduces to the Legendre relation when A = fi = v = and x = r 2 . 

Another generalization of the Legendre relation was given in |AQVV| . 
With the notation 

u = u(r) = F(a — 1, b; c; r), v = v(r) = F(a, b; c; r), 

u\ = u(l — r), v\ = v(l — r), 
the authors considered the function 

L(a, 6, c, r) = uv\ + u\v — vvi, 
proving, in particular, that 

r 2 (c) 



L(a, 1 — a, c, r) 



r(c + a - l)r(c - a + 1) 



This identity reduces to Elliott's formula in case A = v = 1/2 — a and 
\i = c + a — 3/2. In |AQVV| it was conjectured that for a,b £ (0,1), 



a + b < 1(> 1), L(a, b, c, r) is concave (convex) as a function of r on (0, 1). In 
|KaV] Karatsuba and Vuorinen determined, in particular, the exact regions 
of abc— space in which the function L(a, b, c, r) is concave, convex, constant, 
positive, negative, zero, and where it attains its unique extremum. 

In [BNPVJ Balasubramanian, Naik, Ponnusamy, and Vuorinen obtained a 
differentiation formula for an expression involving hypergeometric series that 
implies Elliott's identity. This paper contains a number of other significant 
results, including a proof that Elliott's identity is equivalent to a formula 
of Ramanujan |Bt| p. 87, Entry 30] on the differentiation of quotients of 
hypergeometric functions. 
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5.5. Some approximations for X(r) by arthr. G. D. Anderson, M. K. 
Vamanamurthy, and M. Vuorinen [AVV2j approximated %{r) by the inverse 
hyperbolic tangent function arth, obtaining the inequalities 

,„ „. 7T / arthrX 1 ^ 2 . , ir arthr 

M ~ 2 (—J <X(r)<~ — , 

for < r < 1. H. Alzer and S.-L. Qiu |A1Q2| refined (JSSJ) as 

. „. 7T / arth r \ 3 ^ 4 . . 7r arth r 
(5-7) - < X(r) < -- 



2 V r J 2 r 

with the best exponents 3/4 and 1 for (arth r)/r on the left and right, re- 
spectively. Seeking to improve the exponents in (|5.7j) . they conjectured that 
the double inequality 

7r / arth r \ 3 / 4+Qr yj- / arthr^ 3 / 4+/3r 



(5-8) - j < 0C(r) < - 

holds for all < r < 1, with best constants a = and /3 = 1/4. Very 
recently Y.-M. Chu et al |CWQi| gave a proof for this conjecture. 

S. Andras and A. Baricz |AB| presented some improved lower and upper 
bounds for %{r) involving the Gaussian hypergeometric series. 

5.6. Approximations for £(r). In |GQlj B.-N. Guo and F. Qi have ob- 
tained new approximations for £(r) as well as for %{r). For example, they 
showed that, for < r < 1, 

vr 1 (l + r) 1 ^ „ . . 7T-1 l-r 2 1 + r 



In recent work |CWQJ|rWCl|WCQJ| Y.-M. Chu et al have obtained estimates 
for £(r) in terms of rational functions of the arithmetic, geometric, and root- 
square means, implying new inequalities for the perimeter of an ellipse. 

5.7. Generalized complete elliptic integrals. For < a < min{c, 1} 
and 0<b<c<a + b, define the generalized complete elliptic integrals of 
the first and second kind on [0, 1] by |HVV| 

(5.9) % aAc = 3C 0)6)C (r) = ^|l^F(a, 6; c; r 2 ), 

(5.10) £ aAc = £ aAc (r) = ^^-F(a - 1, 6; c; r 2 ), 
(5-11) 3C' aAc = K aAc (r>) and £^ c = £ a , fe , c (r'), 



for r 6 (0, 1), r' = yl — r 2 . The end values are defined by limits as r tends 
to 0+ and 1 — , respectively. Thus, 

3Ca,6,c(0) = £a,b, C (0) = 

and 

P m l %^(c,c + l-a-i) 

£a,6,c(l = 57 i -i M 1 K o,6,c(l) = OO. 

2 -fc>(c + 1 — a, c — b) 
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Note that the restrictions on the parameters a, b, and c ensure that the 
function % a ^, c is increasing and unbounded, whereas £ a bc i s decreasing 
and bounded, as in the classical case a = b = 1/2, c = 1. 

V. Heikkala, H. Linden, M. K. Vamanamurthy, and M. Vuorinen |HVV| 
IHLVV] derived several differentiation formulas, and obtained sharp mono- 
tonicity and convexity properties for certain combinations of the generalized 
elliptic integrals. They also constructed a conformal mapping sn aj b jC from a 
quadrilateral with internal angles bir, (c — b)n, (1 — a)n, and (1 — c + a)n 
onto the upper half plane. These results generalize the work of |AQVV| . 
For some particular parameter triples (a, b, c) there are very recent results 
by many authors |Bar9l IWZC1I IZWC31 [ZQW]. 

With suitable restrictions on the parameters a, b, c, E. Neuman [N2J has 
obtained bounds for "X a ,b,c an d £ a ,f>,c an d for certain combinations and prod- 
ucts of them. He has also proved that these generalized elliptic integrals are 
logarithmically convex as functions of the first parameter. 

In 2007 A. Baricz [Bar 3] IBar8] established some Turan-type inequalities 
for Gaussian hypergeometric functions and generalized complete elliptic in- 
tegrals. He also studied the generalized convexity of the zero-balanced hy- 
pergeometric functions and generalized complete elliptic integrals [Bar5j (see 
also |Bar9l IBarfl iBlrg] ). 

5.8. The generalized modular function and generalized linear dis- 
tortion function. Let a,b,c > with a + b > c. A generalized modular 
equation of order (or degree) p > is 

F(a,b;c;l-s 2 ) F(a, bj c; 1 - r 2 ) 
F{a,b;c; s z ) F(a,b;c;r z ) 

The generalized modulus is the decreasing homeomorphism fi a ,b,c ■ (0, 1) — > 
(0,oo), defined by 

(5.13) ^IrWI'^l-^l 



2 F(a,b;c;r 2 ) 
The generalized modular equation (|5.12p can be written as 

Va,b,c( S ) = Pl J 'a,b,c( r )- 

With p = 1/K, K > 0, the solution of f|5 . 12[) is then given by 

s = (p a ^ b ' c (r) = ^~l^a,b,c{r)/K). 
Here (f'^ ,,c is called the (a, b,c)- modular function with degree p = 1/K 



|AQVV[ IHW1IHEW] . Clearly the following identities hold: 



^a,fe, C (r)/X a ,6,c(?' / ) 



B(a,b) 



2 



a.b.cr \2 i a,b,c/ i\2 1 

<Pk v) + <Pk v ) = L 
In |HLVV] . the authors generalized the functional inequalities for the modu- 
lar functions and Grotzsch function /U proved in |AQVV| to hold also for the 
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generalized modular functions and generalized modulus in the case b = c — a. 

For instance, for < a < c < 1 and K > 1, the inequalities 

(5.14) 

^ c _ a , c (l - V(l -«)(!-«)) < ^ c ~ aM ± ^ < /io^cCVS) 
hold for all u,t £ (0, 1), with equality if and only if it = i, and 

(5.15) r l / K < <p a /- a ' c (r) < e (i-i/*)«(«.c-«)/2 r l/tf ; 

(5.16) > ^ c " a,c (r) > e (l-J0*(«M~>)/2 r *-. 

For the special case of o = 1/2 and c = 1 the readers are referred to |AVV5| . 
G.-D. Wang et al |WZC1| presented several sharp inequalities for the gen- 
eralized modular functions with specific choice of parameters c = 1 and 
b = 1 - a. 

A linearization for the generalized modular function is also presented in 
|HLVV| as follows: Let p : (0, 1) — > (—00,00) and q : (—00,00) — > (0, 1) be 
given by p(x) = 21og(x/x') and q(x) = p~ l {x) = W e x /(e x + 1), respectively, 
and for a G (0, 1), c G (a, 1], K G (1, 00), let g, /i : (—00, 00) — > (—00, 00) be 
defined by <?(x) = p(ip^f a ' c (q(x))) and /i(x) = p{ip a ^J' K a,c {q{x))). Then 

, , / Ifx, if x > 0, , , , , / x/if, if x > 0, 

o(x) > < . ' , r ' and nix) < < ' .„ _ 

aw - I x/K, if x < 0, I Kx -. lf x < 0- 

In the same paper the authors also studied how these generalized functions 
depend on the parameter c. Corresponding results for the case c = 1 can be 
found in the articles |AQVV[ [QV2| IZWClj . 

Recently B. A. Bhayo and M. Vuorinen |BV2| have studied the Holder 
continuity and submultiplicative properties of (p^l ' c (r) in the case where c = 
1 and b = 1 — a, and have obtained several sharp inequalities for ip^ ~ a ' (r). 

For x,K G (0,oo), define 



VkW = 77 1 s = l Pk 



a,l— a,l/ \ IX 



r 



1 + x' 



*(MQ=1 M-a,i //^ I ■ Ma,l) = l. 



and the generalized linear distortion function 

For a = 1/2, these two functions reduce to the well-known special case de- 
noted by t]k{x) and A(-fT), respectively, which play a crucial role in quasicon- 
formal theory. Several inequalities for these functions have been obtained as 
an application of the monotonicity and convexity of certain combinations of 
these functions and some elementary functions; see |BV21 CWZQ MQZC1 



MQZC2 


ZWUQJ. 


in |CWZ 


Q : for a 



2n a K (x)n a K (y) ^ a ( 2xy 



\vk( x ) + VkW \ x + yJ\ 

f rj a K (x) +Ti a K {y) a f x + y 



< \l if K {x)r] a K {y) <min| ,rf K 



TOPICS IN SPECIAL FUNCTIONS III 



25 



with equality if and only if x = y. 

6. Inequalities for power series 

The following theorem is an interesting tool in simplified proofs for mono- 
tonicity of the quotient of two power series. 

6.1. Theorem. fHWl Theorem 4.3] Let £~ =0 a n x n and^ =Q b n x n be two 
real power series converging on the interval (-R, R). If the sequence {a n /b n } 
is increasing ( decreasing ) and b n > for all n, then the function 

Eoo n 
n=0 a n x 



/(*) 



En=0 b n 



X" 



is also increasing (decreasing) on (0,R). In fact, the function 

oo 

f(x)C£b n x n ) 2 

n=0 

has positive Maclaurin coefficients. 

A more general version of this theorem appears in |BK| and [FY] Lemma 
2.1]. This kind of rule also holds for the quotient of two polynomials instead 
of two power series. 

6.2. Theorem. |HVVI Theorem 4.4] Let f n (x) = Y2=o a kX k and g n (x) = 
s ^Jk = Qbj i x k be two real polynomials, with b^ > for all k. If the sequence 
{o-k/bk} is increasing (decreasing), then so is the function f n (x)/g n (x) for 
all x > 0. In fact, g n f' n ~ fn9n ^ as positive (negative) coefficients. 

In 1997 S. Ponnusamy and M. Vuorinen |PV] refined Ramanujan's work 
on asymptotic behavior of the hypergeometric function and also obtained 
many inequalities for the hypergeometric function by making use of Theorem 
16.11 Many well-known results of monotonicity and inequalities for complete 
elliptic integrals have been extended to the generalized elliptic integrals in 

[HWllHEY]- 

Motivated by an open problem of G. D. Anderson et al |AVV3j . in 2006 
A. Baricz |Bar2| considered ratios of general power series and obtained the 
following theorem. Note the similarity with the inequality (|5.14[) . 

6.3. Theorem. Suppose that the power series f(x) = Yln=o a n xTl with a n > 
for all n > is convergent for all x £ (0,1), and also that the sequence 
{(n + l)a n +i/a n — n} n >o is strictly decreasing. Let the function rrif : (0, 1) — > 
(0, oo) be defined as mj-(r) = /(l — r 2 )/f(r 2 ). Then 



^Y[mf(ri) <m f ■ ^ 




for all ?r,r2, • • • , rjt S (0, 1), where equality holds if and only if r± = r<i 
■ ■ ■ = rfc. In particular, for k = 2 the inequalities 



^mf(ri)mf(r 2 ) < mj( A /rir 2 ) 



+ > 



m/(n) m f (r 2 ) m f (^/ryr^)' 
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m/(ri) + m/(r 2 ) > 2m/ ^1 - - rf)(l - r~f)\ 

hold for all T\,Ti 6 (0>1); aric ^ W a ^ i/iese inequalities equality holds if and 
only if n = T2 . 

The following Landen-type inequality for power series is also due to Baricz 
[BarT] . 

6.4. Theorem. Suppose that the power series f(x) = EnLo a nX n with a n > 
for all n > is convergent for all x € (0, 1), and i/iai for a given 5 > 1 i/ie 
sequence {n!a n /(log <5) n } n >o is decreasing. If A/(x) = f{x 2 ), then 

A ' (r£) < ^ (r) 

/io/ds /or a// r € (0, 1) and p > o" 4 ^ 5 . 

G. D. Anderson, M. K. Vamanamurthy, and M. Vuorinen |AVV9| studied 
generalized convexity and gave sufficient conditions for generalized convexity 
of functions defined by Maclaurin series. These results yield a class of new 
inequalities for power series which improve some earlier results obtained by 
A. Baricz. More inequalities for power series can be found in [Bar9, CWZQ|. 



In 1928 T. Kaluza gave a criterion for the signs of the power series of a 
function that is the reciprocal of another power series. 

6.5. Theorem. |K] Let f(x) = Xm>o a nX n be a convergent Maclaurin series 
with radius of convergence r > 0. If a n > for all n > and the sequence 
{O"n}n>o is log-convex, that is, for all n> 

(6.6) a 2 n < a n -ia n+ i, 

then the coefficients b n of the reciprocal power series 1/ f(x) = Yln>o^"-X n 
have the following properties: bo = 1/eto > and b n < for all n > 1. 

In 2011 A. Baricz, J. Vesti, and M. Vuorinen [BVVJ showed that the 
condition (|6.6p cannot be replaced by the condition 



a n-l + a n+l\ ^ 



2 J 1 

for any t > 0. However, it is not known whether the condition (|6.6p is 
necessary. 

In 2009 S. Koumandos and H. L. Pedersen [KPJ proved the following in- 
teresting result, which deals with the monotonicity properties of the quotient 
of two series of functions. 

6.7. Theorem. |KP| Lemma 2.2] Suppose that at > 0, bk > and that 
{uk(x)} is a sequence of positive C 1 -functions such that the series 

oo oo 

^ afcttfc (x) and ^huf' (z), 1 = 0,1, 
k=o k=0 
converge absolutely and uniformly over compact subsets of [0, oo). Define 

Efclo a kUk{x) 



J2T=o b ku k (x)' 
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(1) If the logarithmic derivatives u' k (x)/uk{x) form an increasing sequence of 
functions and if a^/bk decreases (resp. increases) then f(x) decreases (resp. 
increases) for x > 0. 

(2) If the logarithmic derivatives u' k (x)/uk(x) form a decreasing sequence of 
functions and if ak/bk decreases (resp. increases) then f(x) increases (resp. 
decreases) for x > 0. 



7. Means 

A homogeneous bivariate mean is defined as a continuous function M : 
M + x M + — > M satisfying min{x, y} < M,(x, y) < max{x, y} and M(Ax, Ay) = 
AM(x, y) for all x, y, A > 0. Important examples are the arithmetic mean 
A(a,b), the geometric mean G(a,b), the logarithmic mean L(a,b), the iden- 
tric mean I(a, b), the root square mean Q(a, 6), and the power mean M r (a, b) 
of order r, defined, respectively, by 

b) = G(a, b) = Vol, 



1 / a -\ V(^) 



L(a, b) = — -, I(a, b) = - 7r 

v ' ; loga-log6' v ' ; e 1 " 



2 -I- h 2 r n r b r 

Q(a,b) = d^—, M r (a,b) = 



7.1. Power means. The weighted power means are defined by 
M x (io;a,b) = (ua x + (1 - w)& A ) 1/A (A / 0), 

Mq(uj; a, b) = a^b 1 ^^, with weights u, 1 — ui > 0. The power means are the 
equally-weighted means M\(a, b) = M\(l/2; a, b). As a special case, we have 
M (l/2;a,b) = G(a,b). 

In [ K2] O. Kouba studied the ratio of differences of power means 

Mj'(a,b)-GP(a,b) 

KM,p;a,&) = MfM) _ GPM) , 

finding sharp bounds for p(s,t,p; a,b) in various regions of sip-space with 
a, 6 positive and a ^ b. This work extends results of H. Alzer and S.-L. 
Qiu |AIQT1 , T. Trif (Tf], O. Kouba (ED, S.-H. Wu |Wu], and S.-H. Wu and 
L. Debnath [WD3J. O. Kouba also extended the range of validity of the 
following inequality, due to S.-H. Wu and L. Debnath [WD3J: 

2-p/r _ 2-p/ s (a, b) - Mj (a, b) r-s 

2-p/t _ 2-p/» < Mf(a,b) - Mi(a,b) < t-s 

to the set of real numbers r,t,s,p satisfying the conditions < s < t < r 
and0<p< (4t + 2s)/3. 
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7.2. Toader means. If p : R + — > M + is a strictly monotonic function, then 
define 

1 

2^ 



p{(a n cos 2 6 + b n sin 2 Of^alO if n / 0, 
— f n 2 *" p(a cos2 e 6 sin2 if n = 0, 

*« 27T 



f(a,b;p,n) 



where a, 6 are positive real numbers. The Toader mean [ToJ of a and b is 
defined as T(a,b;p,n) = p~ 1 (f(a,b;p,n)). It is easy to see that the Toader 
mean is symmetric. For special choices of p, let Tg >ri (a, b) = T(a,b;p,n) 
if p(x) = x q with 9 / 0, and To >n (a, 6) = T(a,b;p,n) if = logx. 

The means Tq n belong to a large family of means called the hypergeometric 
means, which have been studied by B. C. Carlson and others [BC, IC1| ICT] , 
It is easy to see that T q)Tl is homogeneous. In particular, we have 

To t2 (a,b) = A(a,b), 

T- 2j2 (a,b) = G(a,b), 

T 2 , 2 {a,b) = Q{a,b). 

Furthermore, the Toader means are related to the complete elliptic integrals: 
for a > b > 0, 

TLi 2 (a,o) 



2^(^l-(fe/a) 2 ) 
and 

2o 



T h2 (a,b) = — £(^l-(6/a) 2 ). 

7T 



In 1997 S.-L. Qiu and J.-M. Shen |QS| proved that, for all a,b > with 
a ^ b, 

M 3/2 (a,b) <T 1)2 (a,6). 

This inequality had been conjectured by M. Vuorinen |V]. H. Alzer and S.-L. 
Qiu |A1Q1| proved the following best possible power mean upper bound: 

Ti )2 (a,b) < M log 2/ io g (7r/2) (a,b). 

Very recently, Y.-M. Chu and his collaborators fCWTl ICW21 |CWQ| obtained 
several bounds for T\ 2 with respect to some combinations of various means. 

7.3. Seiffert means. The Seiffert means S\ and £2 are defined by 

S l( a > b ) = n a ■ b n -h > a ^ S l( a ' a ) = a ' 
2 arcsm ^-j 

a+b 

and 

5*2(0, 6) = r, a ^ b, S 2 (a, a) = a. 

2arctan^ 

It is well known that 

mi<L< 2G+A 



3 

Sandor proved similar results for Seiffert means |S1| IS2] : 

G + 2A 



(7.1) ^G <Sl< ^-^<I 
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and 

(7.2) 3/Q2A <S2 <^9.. 

The inequalities (|7.ip and (|7.2p are special cases of more general results 
obtained by E. Neuman and J. Sandor |NS1| [NS2J. 

7.4. Extended means. Let a, b 6 (0,oo) be distinct and s,t G M\ {0}, 
s 7^ t. We define the extended mean [St] with parameters s and t by 



E s , t (a,b) 
and also 



6* 



, . f a s log a - 6 s log 6 
£ s , s (a, 6) = exp ( - + — 

1/8 



E s0 (a,b)=[— — ^—r-r) and E 00 (a,b) = Vab. 

We see that all the classical means belong to the family of extended means. 
For example, E21 = A, Eq } q = G, £Li _2 = H, and = L, and, more 
generally, M\ = £^2A,A f° r A G M. The reader is referred to the survey |Q] for 
many interesting results on the extended mean. 

In 2002 P. A. Hasto [HalJ studied a certain monotonicity property of ratios 
of extended means and Seiffert means, which he called a strong inequality. 
These strong inequalities were shown to be related to the so-called relative 
metric |Ha2| Ha3j. 

7.5. Means and the circular and hyperbolic functions. It is easy to 
check the following identities: 

(7.3) A(l + sinx, 1 — sinx) = 1, G(l + sinx, 1 — sinx) = cosx, 



(7.4) Q(l -l-sinx, 1 — sinx) = y 1 + sin 2 x, 5i(l + sinx, 1 — sinx) 



smx 



.r 



(7.5) A(e x ,e~ x ) = coshx, G(e x ,e~ x ) = l, Q(e x ,e~ x ) = Vcosh2x, 

(7.6) L(e x ,e- X ) = ^1 I( e x ,e~ x ) = e xcothx ~ l , 

x 

(7.7) S 1 (e x ,e- X ) = Si " h \ . , S 2 {e x ^ Huh ''' 



arcsin(tanhx) ' ' arctan(tanhx) 

One can get many inequalities by combining the above identities and in- 
equalities between means. For example, combining (|7.ip and (|7.4p we have 

sin x cos x + 2 



v^cosx < < 

x 6 

where the second inequality is the well-known Cusa-Huygens inequality, and 
combining (|7.2p . (|7.5p . and (|7.7p . we have 

sinh x cosh x + 2\/cosh 2x 



v/ (cosh 2x)(cosh x) < r < 

arctan(tanhx) 3 

More inequalities on mean values and trigonometric and hyperbolic functions 
can be found in [S4, N9, YalJ and references therein. 
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7.6. Means and hypergeometric functions. In 2005 K. C. Richards jRi] 
obtained sharp power mean bounds for the hypergeometric function: Let 
< a, b < 1 and c > max{— a, b}. If c > max{l — 2a, 2b}, then 

M A (1 - b/c; 1, 1 - r) < F(-a, b; c; r) 1/a 

if and only if A < f±f. If c < min{l - 2a, 2b}, then 

M M (1 - b/c; 1, 1 - r) < F{-a, b; c; rfl a 

if and only if /i > y4^. These inequalities generalize earlier results proved 
by Carlson |C2] . 

For hypergeometric functions of form F(l/2 — s, 1/2 + s; 1; 1 — r^) 9 , J. M. 
Borwein et al |BBG| exhibited explicitly iterations similar to the arithmetic- 
geometric mean. R. W. Barnard et al [BRJ presented sharp bounds for 
hypergeometric analogs of the arithmetic-geometric mean as follows: For 
< a < 1/2 and p > 0, 

M x {a; 1, r) < F(a, 1 - a; 1; 1 - r p )" 1/(ap) < M^a; 1, r) 

if and only if A < and ^ > p(l — a)/2. 

Some other inequalities involving hypergeometric functions and bivariate 
means can be found in the very recent survey [BRTJ. 

For any two power means M\ and M^, a function / is called M\^- convex 
if it satisfies 

f(M x (x,y))<M ll (f(x),f(y)). 

Recently many authors have proved that the zero-balanced Gaussian hy- 
pergeometric function is M^A-convex when A € { — 1,0,1}. For details see 
[BPV] . |AVV9] . [Ba79] . and |CWZQ| . A. Baricz |Bar5] generalized these re- 
sults to the M^A-convexity of zero-balanced Gaussian hypergeometric func- 
tions with respect to a power mean for A G [0, 1]. X.-H. Zhang et al |ZWC2| 
extended these results to the case of MA jA1 -convexity with respect to two 
power means: For all a, b > 0, AG (— oo, 1], and fj, G [0, oo) the hypergeo- 
metric function F(a, b;a + b; r) is MA^-convex on (0, 1). 

The following interesting open problem is presented by A. Baricz [Bar8j: 

Open problem. If mi andrri2 are bivariate means, then find conditions on 
o-i, «2 > and c > for which the inequality 

mi {F ai {r),F a2 {r)) < (>) 

-^m 2 (ai,a 2 ) ( r ) 

holds true for all r G (0, 1), where F a (r) = F(a, c — a;c;r). 

7.7. Means and quasiconformal analysis. Special functions have always 
played an important role in the distortion theory of quasiconformal map- 
pings. G. D. Anderson, M. K. Vamanamurthy, and M. Vuorinen |AVV5| 
have systematically investigated classical special functions and their exten- 
sive applications in the theory of conformal invariants and quasiconformal 
mappings. Some functional inequalities for special functions in quasicon- 
formal mapping theory involve the arithmetic mean, geometric mean, or 
harmonic mean. For example, for the well-known Grotzsch ring function 
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/x and the Hersch-Pfluger distortion function tpx, the following inequalities 
hold for all s, t G (0, 1) with s ^ t: 



and 

■sj ipK(s)y>K(t) < ipK(Vsi) for K > 1. 
Recently, G.-D. Wang, X.-H. Zhang, and Y.-M. Chu |WZC21 IWZC3| have 
extended these inequalities as follows: 

M x (n(s),n(t)) < fj,(M x (s,t)) if and only if A < 0, 

M\(ipK(s),(pi<{t)) < <pK(M\(s,t)) if and only if A > and K > 1, 

and 

M x (tpK(s),(p K (t)) > tp K (M x (s, t)) if and only if A > and < K < 1. 

Some similar results for the generalized Grotzsch function, generalized 
modular function, and other special functions related to quasiconformal anal- 
ysis can be found in |QQWG[ |WQZG[ IWZG41 IWZJ11 IWZJ2] . 



8. Epilogue and a view toward the future 

In earlier work we have listed many open problems. See especially |AVV1| 
pp. 128-131] and [AVV5, p. 478]. Many of these problems are still open. In 
Sections 4, 6, and 7 above we have also mentioned some open problems. 

Finally, we wish to suggest some ideas for further research. In a frequently 
cited paper |Lt| P. Lindqvist introduced in 1995 the notion of generalized 
trigonometric functions such as sin p , and presently there is a large body 
of literature about this topic. For the case p = 2 the classical functions are 
obtained. In 2010 R. J. Biezuner et al [BEM| developed a practical numerical 
method for computing values of sin p . Recently, S. Takeuchi [Taj has gone 
a step further, introducing functions depending on two parameters p and q 
that reduce to the p-functions of Lindqvist when p = q. In |BV1| IBV3] the 
authors have continued the study of this family of generalized functions, and 
have suggested that many properties of classical functions have a counterpart 
in this more general setting. It would be natural to generalize the properties 
of trigonometric functions cited in this survey to the (p, g)-trigonometric 
functions of Takeuchi. 
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